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ABSTRACT 



(-1 

• ^ ^1 We investigate the structural, kinematical, and chemical properties of globular 

I . cluster systems (GCSs) in galaxies of different Hubble types in a self-consistent man- 

O ' ner based on high-resolution cosmological N-body simulations combined with semi- 

analytic models of galaxy and globular cluster (GC) formation. We focus on corre- 
^ i lations between the physical properties of GCSs and those of their host galaxies for 

' ^ 10^ simulated galaxies located at the centres of dark matter halos (i.e. we do not 

consider satellite galaxies in sub- halos). Our principal results, which can be tested 
against observations, are as follows. The majority (~ 90%) of GCs currently in halos 
^ ' are formed in low-mass galaxies at redshifts greater than 3 with mean formation red- 

CN ! shifts of z = 5.7 (12.7 Gyrs ago) and 4.3 (12.3 Gyrs ago) for metal-poor GCs (MFCs) 

' and metal-rich GCs (MRCs) respectively. About 52 % of galaxies with GCs show clear 

bimodality in their metallicity distribution functions, though less luminous galaxies 
with Mb fainter than —17 are much less likely to show bimodality owing to little or no 
' MRCs. The number fraction of MRCs does not depend on Hubble type but is generally 

I smaller for less luminous galaxies. The specific frequencies (^n) of GCSs are typically 

higher in ellipticals (S'n ~ 4.0) than in spirals (S'n ~ 1-8), and higher again (S'n ~ 5.0) 
for galaxies located at the centers of clusters of galaxies. The total number of GCs per 
unit halo mass does not depend strongly on Mb or Hubble type of the host galaxy. 
The mean metallicities of MFCs and MRCs depend on Mb such that they are higher 
. ^ in more luminous galaxies, though the dependence is significantly weaker for MFCs. 

^ . The spatial distributions of MRCs are more compact than those of MFCs and we find 

that the half-number radii of MFCs (re,mpc) correlate with the halo masses (Mh) such 
that r-cmpc oc Afh"'^*- There is no significant difference in velocity dispersions between 
MFCs and MRCs. We qualitatively compare our results to observational data where 
possible. Finally, we discuss these results in the wider context of galaxy formation and 
evolution. 

Key words: globular clusters: general - galaxies: star clusters - galaxies: evolution - 
galaxies:stellar content 
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1 INTRODUCTION 

The physical properties of GCSs in galaxies have long been 
considered to be "fossil records" that contain vital infor- 
mation on galaxy formation and evolution (e.g., Searle & 
Zinn 1978; Harris 1991; Ashman & Zepf 1998; West et al. 
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2004). For example, the observed bimodal colour distribu- 
tions, higher specific frequencies (S'n) and higher fraction 
of metal-poor GCs (MPCs) in elliptical galaxies has moti- 
vated various formation scenarios for elliptical galaxies (e.g.. 
Ashman & Zepf 1992; Forbes et al. 1997; Cote et al. 1998). 
Correlations between physical properties of GCSs and those 
of their host galaxies also have provided clues to the bet- 
ter understanding galaxy formation For example, Strader 
et al. (2004) showed that the mean colours of MPCs and 
MRCs correlate with the luminosity (L) of their host galaxy. 
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From a survey of Virgo cluster carly-typo galaxies, Peng et 
al. (2006, P06) calculated from g-z colours, the relations 
[Fe/HJg^ oc lO i6±o.o4 ^PCs and [Fe/U]^^ oc L°-^«±° °2 
for MRCs. 

In addition to the above chemical property of GCSs, the 
kinematical and structural properties of GCSs have been in- 
vestigated both for MFCs and MRCs (Kissler-Patig & Geb- 
herdt 1998; Cote et al. 2001; Rhode & Zepf 2004; Richtler 
et al. 2004 Peng et al. 2005; Bergond et al. 2006; Bridges et 
al. 2006; Pierce et al. 2006; Romanowsky 2006; Hwang et al. 
2007; Woodley et al. 2007). Feng et al. (2005) showed that 
the GCS in NGC 5128 has a significant amount of global 
rotation whereas Richtler et al. (2004) did not find any sig- 
nificant rotation in the GCS of NGC 1399, which suggests 
a great diversity in the kinematics of GCSs in early-type 
galaxies. The slopes of power-law density profiles of GCSs 
are also observed to be diverse and only loosely correlated 
with the luminosity of their host galaxy (e.g., Harris 1986; 
Ashman & Zepf 1998). 

In comparison with this remarkable progress in observa- 
tional studies of GCSs (Brodie & Strader 2006), theoretical 
modeling has not been as well developed as to provide useful 
predictions that can be compared with the wide range of ob- 
servations on GCS properties. This is mainly because, both 
subpc-scale formation processes of GCs in galaxies and kpc- 
and Mpc-scales merging/interaction processes over a Hubble 
time need to be modeled in a fully self-consistent manner in 
order that the GCS projjcrties and GCS-host relations can 
be examined quantitatively. 

Beasley et al. (2002; B02) first investigated the phys- 
ical properties of GCSs and their correlations with those 
of host galaxies based on a semi-analytic model of hierar- 
chical galaxy formation in a ACDM universe. B02 however 
could not discuss the structural and kinematical proper- 
ties of GCSs owing to the limitations of the adopted semi- 
analytic model. Although previous collisionlcss N-body sim- 
ulations with GCs have provided some reasonable explana- 
tions for the origin of structural and kinematical properties 
of GCSs (e.g., Bekki et al. 2005; Bekki & Forbes 2006), they 
did not model the chemical properties of GCSs. Kravtsov & 
Gnedin (2005) carried out a high resolution gas + N-body 
simulation of GC formation in a ACDM universe. In their 
model of a Milky Way like galaxy, GCs first formed at red- 
shift ~ 12 with peak formation occurring at z ~ 4. However 
they only simulated MFCs and only from formation to a 
rcdshift of 3. Thus no sirrmlations have yet addressed the 
structural, kinematical, and chemical properties of GCSs in 
a self-consistent cosmologically motivated simulation from 
formation to the current epoch. 

The purpose of this paper is to investigate both the dy- 
namical and chemical properties of GCSs based on cosmolog- 
ical N-body simulations with semi-analytic models of galaxy 
formation and thereby compare the results with the above- 
mentioned growing number of observational constraints on 
GCS-host relations. We focus particularly on the follow- 
ing physical properties of GCSs in galaxies : (i) GC metal- 
licity distribution functions (MDFs), (ii) number fractions 
of MRCs (/mrc), (iii) number fractions of GCS with bi- 
modal MDFs (/bimo), (iv) specific frequencies of GCs (Sn), 
(v) mean GCS metallicities ([Fe/H]), (vi) GCS half-number 
radii (re), and (vii) GCS velocity dispersions (a). We mainly 
investigate correlations between these GCS properties and 



their host properties such as Mb and Hubble type. Here 
wo do not discuss other important observational results of 
GCSs, such as the blue tilt (e.g., Strader et al. 2006) or the 
GC luminosity function dependence on galaxy luminosities 
(e.g., Jordan et al. 2006). Some of these aspects have been 
already discussed in our previous papers (e.g., Bekki et al. 
2007a, B07). Here we also restrict our analysis to central 
galaxies and do not consider satellites in sub-halos. 

The plan of this paper is as follows: In the next section 
wc describe our rmmerical method of cosmological N-body 
simulations, semi-analytic model of galaxy formation, and 
formation model for GCs. In §3, we present our numerical 
results on GCS properties in galaxies with different lumi- 
nosities and Hubble types. In §4, we discuss our results of 
GCS-host relations in the context of galaxy formation and 
evolution. We summarize our conclusions in §5. 



2 THE MODEL 

We consider that GCs can be formed within any galaxy at 
any redshift, if the physical conditions required for GC for- 
mation are satisfied in the galaxies. We therefore investi- 
gate GC formation rates in any virialized halo in a high- 
resolution cosmological simulation based on a ACDM cos- 
mology model. The physical properties of GCs (e.g., metal- 
licities) are determined by their host galaxies at the epoch 
of their formation. GCs formed in low-mass galaxies at high 
z can be tidally stripped during the hierarchical merging 
of galaxies to finally become GCs within a giant galaxy 
at « = 0. Since numerical methods and techniques of the 
present cosmological simulations and those of semi-analytic 
models of galaxy formation have been given in our previ- 
ous papers (e.g., Yahagi et al. 2004; Nagashima et al. 2005, 
N05), we only briefiy describe them in the present study. 
We here focus on (i) how to identify GCs in virialized dark 
matter halos and (ii) how to allocate physical properties to 
the identified GCs based on their host galaxy properties. 

2.1 Simulations 

We simulate the large scale structure of GCs in a ACDM 
Universe with fi = 0.3, A = 0.7, Ho = 70 km s"^ Mpc"\ 
and as = 0.9 by using the Adaptive Mesh Refinement 
A?"— body code developed by Yahagi (2005) and Yahagi et 
al. (2004), which is a vectorized and parallelized version of 
the code described in Yahagi & Yoshii (2001). We use 512^ 
coUisionless dark matter (DM) particles in a simulation with 
the box size (-Rt) of 70/i~^Mpc and the total mass (A^t) of 
4.08 X 10^'' M©. We start sirrmlations at Zi = 41 and follow 
it until 2: = in order to investigate the physical properties 
of GCs outside and inside of virialized dark matter halos. 
We use COSMICS (Cosmological Initial Conditions and Mi- 
crowave Anisotropy Codes), which is a package of Fortran 
programs for generating Gaussian random initial conditions 
for non-linear structure formation simulations (Bertschinger 
1995, 2001). 

Our method of identifying GCs (or "GC particles") and 
following their evolution is described as follows. Firstly, we 
select virialized dark matter subhalos at a given redshift 
by using the friends-of-friends (FoF) algorithm (Davis et al. 
1985) with a fixed linking length of 0.2 times the mean DM 
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Table 1. Model parameters in the simulation. 

n A //o (km s-i Mpc-i) 0-8 A/t^ ^ ^trun^ 13 f_ SNe feedback _9 

0.3 0.7 70 0.9 4.0" 70~ 4l" 6 ~ O.o" O.oT strong ~ 



^ The total mass of a simulation in units of 10^* Mq . 
^ The box size of a simulation in units of Mpc 
The redshift (2) at which a simulation starts. 

The redshift of truncation of GC formation: GC formation is completely truncated in halos that are virialized later than the redshift. 
^ The coefficient in the functional form describing the dependence of GC formation rate on the mass ratios of merging two galaxies. 
f The coefficient in the functional form describing the dependence of GC formation rate on the gas mass fractions of galaxies. 
9 The details of the strong and weak feedback effects are given in N05. 



particle separation. The minimum particle number A^min for 
halos is set to be 10. For each individual virialized subhalo, 
the central particle is labeled as a "GC" particle. This pro- 
cedure for defining GC particles is based on the assumption 
that energy dissipation via radiative cooling allows baryons 
to fall into the deepest potential well of dark matter halos 
and finally to be converted into GCs and stars. 

The adopted initial distributions of GCs with respect to 
their host halos would be oversimplified, given that possible 
candidates of forming GCs are not necessarily in the central 
regions of galaxies at z = 0. We think that the adopted as- 
sumption can be regarded as reasonable, because previous 
observations suggested that a significant fraction (or even 
all) of the Galactic GCs originate from nuclei of the Galactic 
building blocks at high 2 (e.g., Zinnecker et al. 1988; Free- 
man 1993). We stress that the predicted spatial distributions 
of GCSs in galaxies at 2 = might well weakly depend on 
the adopted initial distributions of GCs within halos, though 
previous simulations suggested that distributions of GCSs 
at 2; = do not depend strongly on the adopted range of 
reasonable initial GC distributions (Yahagi & Bekki 2005). 

Secondly, we follow GC particles until z = and thereby 
derive their locations {x,y,z) and velocities (wx,Vy,i'z). We 
then identify virialized halos at 2; = with the FoF algo- 
rithm and investigate whether each GC is within the virial 
radius (rvir) of a halo. If GCs are found to be within a halo, 
the physical properties of the GCS are investigated. If a 
GC is not in any halo, it is regarded as an intergalactic 
GC. We don't discuss these intergalactic (i.e., intra-group 
and intra-cluster) GCs further in this current paper (see 
Yahagi & Bekki 2005; Bekki & Yahagi 2006 for details of 
such GCs), though a growing number of observations have 
revealed physical properties of these intergalactic GCs (e.g., 
Bassino et al. 2006; Jones et al. 2006). 

Recent theoretical works have suggested that heating 
and gas loss resulting from reionization can severely sup- 
press star formation in low-mass galaxies during reioniza- 
tion (e.g., Susa & Umemura 2004). It is therefore highly 
likely that GC formation is truncated in low-mass galaxies 
that are virialized after the completion of reionization. In 
order to include the effects of the suppression of GC forma- 
tion via reionization on the final properties of the simulated 
GCSs, we adopt the following somewhat idealized assump- 
tion: If a galaxy is virialized after the completion of reion- 
ization (zreion), then GC formation is totally suppressed in 
the galaxy. Therefore, GC particles formed in galaxies with 
«vir < Zreion are not Considered in the physical properties of 
the simulated GCSs. We define the truncation epoch of GC 



formation as ^trun (rather than z^aion for convenience) in the 
present study. Recent quasar absorption-line studies give a 
lower limit of 6.4 for Zreion (Fan et al. 2003). Guided by these 
observations, we investigate the model with Ztmn = 6. 

Thus the GCS of a galaxy is a collection of GCs that 
are formed within low-mass galaxies (i.e., galaxy building 
blocks) embedded in massive dark matter halos virialized 
at high redshifts. The physical properties of the GCS of a 
galaxy therefore depends on star formation histories, chem- 
ical evolution, and merging histories of the building blocks. 
Chemical evolution and star formation histories of building 
blocks for galaxies are derived from the semi-analytic model 
(N05) which are based on the merging histories of DM halos 
derived from our N-body simulation. 



2.2 Semi-analytic model 

The physical properties of a newly formed GC in a viri- 
alized halo at any redshift are determined by those of its 
host galaxy at that redshift. N05 constructed the Numeri- 
cal Galaxy Catalog (i/GC) based on a semi-analytic model 
combined with high-resolution A'^— body simulations. In the 
present study, we adopt the same semi-analytic model as 
that used by N05 which includes various physical processes 
associated with galaxy formation, such as galaxy merging, 
radiative gas cooling, star formation, supernovae feedback, 
and extinction by internal dust. Since the methods and tech- 
niques are given in detail by N05, we do not describe them 
in the present study. 

The semi-analytic model by N05 can reproduce many 
observations reasonably well, such as cold gas mass-to-stellar 
luminosity ratios of spiral galaxies, faint galaxy number 
counts, cosmic star formation rates, the TuUy-Fisher relation 
for bright spiral galaxies, luminosity functions of local galax- 
ies, and colour-magnitude relations for massive and dwarf el- 
liptical galaxies. The successes and limitations of this model 
in explaining observed galactic properties are given in N05. 

Here we present the results for the properties of GCs 
from our semi-analytic model and compare them where pos- 
sible to observations. We focus on results from the model 
with the "strong feedback effects" (N05) in which galaxies 
with smaller rotational velocities are much more strongly in- 
fluenced by heating. An important caveat in this work is that 
only galaxies located at the centers of halos are analyzed in 
the present study, i.e. we do not investigate satellite galax- 
ies. Thus we effectively focus on brightest cluster, brightest 
group and isolated galaxies in this current work. Observa- 
tional studies are often dominated by non-central galaxies 
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Figure 1. The distribution of GCs projected onto tlie x-y plane 
at z = 0. Here only GC particles that are at the very centers 
of GCSs for galaxies at 2: = are shown for convenience. Thus 
the distribution describes the large-scale distribution of GCSs of 
galaxies in the universe at z = 0. 

(satellites in this context) and this should be born in mind 
when we compare our model predictions to the observations 
below. 

2.3 Formation efficiencies of GCs 

We convert the star formation rate (SFR) in a galaxy into a 
GC formation rate (GCFR) based on the physical properties 

of the host galaxy at a given rcdshift. The vast majority of 
stars arc observed to form in star clusters (SCs) embedded 
within GMCs (Lada & Lada 2003). Strongly bound SCs can 
evolve into GCs or old open clusters whereas weakly bound, 
low-mass ones can be disintegrated into field stars of the host 
galaxy. The mass fraction of new SCs that evolve into GCs 
relative to all new SCs is denoted as CeS and is assumed to 
be dependent on the physical properties of their host galaxy. 
Thus we define the GCFR as: 

GCFR = Ceff X SFR, (1) 

where SFR is an output from the semi-analytic model. 

We consider both observational results by Larsen & 

Richtler (2000, LROO) and simulation ones by Bekki ot al. 
(2002, BFBC02) in order to determine Ccff in a physically 
reasonably way. LROO investigated correlations between spe- 
cific i7— band cluster luminosities Ti^{U) of GCSs and their 
host galaxy properties both for apparently isolated galaxies 
and for interacting/merging ones. LROO found that (i) Ti^iU) 
correlates with SFR per unit area, stellar surface brightness, 
and HI surface density, and (ii) Tl{U) is more than an or- 
der of magnitude higher in strongly starbursting mergers 
(e.g., NGC 1705) than apparently normal galaxies that are 
forming young SCs (see Tabic 1 in LROO). These results by 
LROO imply that CeS can be higher in gala^xies with higher 
surface densities and higher gas mass fractions for a given 
galaxy luminosity. 




redshift 



Figure 2. Cosmic evolution of GC formation rates (GCFRs) as 
a function of redshift (2) for MFCs (solid) and MRCs (dotted). 
The GCFRs are averaged over the volume of the simulation. 




redshift 

Figure 3. Total numbers of GCs formed at redshift z for MFCs 
(solid) and MRCs (dotted). The numbers are normalized by the 
maximum values for < z < 12. 

Numerical simulations of GC formation in merging 
galaxies (BFBC02) found that Ceff depends on the mass 
ratios of the merging spirals (/m or 7712) and their gas mass 
fractions (/g) in such a way that Cefi is higher in mergers 
with larger m2 and /g. These results by BFBC02 combined 
with LROO imply that Ceff is a minimum for isolated galaxies 
with m2 = and maximum for equal mass galaxy mergers 
with m2 — 1. BFBC02 also found that Cog is lower in low 
surface brightness galaxies for a given galaxy mass. Galax- 
ies with their halos virializcd at higher z have more compact 
disks and thus higher stellar densities for a given mass, bary- 
onic mass fraction, spin parameter, and halo circular veloc- 
ity in galaxy formation models based on ACDM (e.g.. Mo et 
al. 1998). The above result by BFBC02 combined with that 
by Mo et al. (1998) therefore strongly suggests that CcB is 
likely to be higher for galaxies formed at higher z. Thus Ces 
is a function of the gas fraction, the merger mass ratio and 
formation rcdshift. 

Wc estimate Cch as follows: 

Ceff- = CoF,(fg)F,Mm)F.{z). (2) 

Wc need to choose the forms of tiicsc tiircc functions 
(Fg, Fm, and Fz) so that they are consistent with previous 
results by LROO and BFBC02. The normalization factor Co 
is determined as follows: 

Co — Fg^max Fni,max Fz^max , (3) 

where Fg^max, Fm.max, and Fz^max are maximum values of the 
above three functions. This method of normalization ensures 
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that Ceff is always equal to or less than 1. If we adopt a 
value of Co significantly smaller than the one given in the 
above equation, physical properties of the simulated GCSs 
(e.g., total numbers of GCs in galaxies) can be much less 
consistent with observations. 

In the present study, we choose elementary functions 
for Ceff. Considering that (i) the dependence of Ces on /m 
appears to be non-linear (BFBC02) and (ii) Ces is much 
higher in major mergers (LR02 and BFBC02), we assume 
the following: 



Fmifm) = !-(! + a)-' + (1 + a - Uy\ 



(4) 



where q is a parameter. This formula ensures that -F'm(/m) 
has a minimum 1) at /m = corresponding to isolated 
disk galaxies and a maximum at /m = 1 corresponding to 
major mergers. Furthermore, if a small value of a (< 0.1) 
is adopted, Fm is more than 10, which corresponds to a 
strongly starbursting major merger in LROO. We adopt the 
same functional form for _Fg(/g): 



^'g(/g) = l-(l+/3)"' + (l + /3-/g 



Thus, a and f3 are free parameters which are chosen so that 
observations of GCSs can be solf-consistontly explained. 

Mo et al. (1998) showed that the sizes (Ra) of disk 
galaxies formed at redshift z are inversely proportional to 
H{z)/Ho for given disk and halo properties (e.g., disk mass 
fraction and halo circular velocity). Therefore Rd can be 
more compact for disk galaxies formed at higher redshifts 
so that mean stellar densities Eg within R^ (Eg oc Rd~^ for 
a given stellar disk mass) are higher. Considering the above 
dependences of disk properties on H(z)/Ho (Mo et al. 1998), 
we adopt the following functional form: 



F,{z) = {H{z)/Hof 



(6) 



We estimate -Fz(2) for each halo at a given redshift based on 
the virialization redshift (zvir) of the halo. 

Here we adopt a simpler yet qualitatively reasonable de- 
pendence of F:i,{z): we do not directly derive the surface mass 
densities by assuming initial sizes of halos, spin parameters, 
and sizes of galaxies based on the results of the SAM by 
N05. Although wo can not discuss the importance of these 
physical properties (e.g., initial spin parameters) in GC for- 
mation in the present model, we can show the importance 
of surface mass (or gas) densities (i.e., the epoch of virial- 
ization) in GC formation of galaxies in a clearer and more 
straightforward way thanks to the adopted model. It should 
be stressed here that without introduction of the F^{z), most 
GCS properties can not be reproduced well in the present 
model. 

Although the adopted functional forms of Fm, Fg, and 
F^ are reasonable at least qualitatively, no observations have 
been carried out which allow us to determine whether the 
adopted functional form is quantitatively consistent with 
observations. Accordingly other functional forms could be 
adopted so that observational properties of GCSs can be re- 
produced s well. Here we do not discuss the results of models 
with difi^crcnt functional forms of Fm, F^, and Fz- 

Tl(U) in some merging galaxies (e.g., NGC 1705) are 
observed to be as large as 15, which is a factor of ~ 20 
larger than the average value (= 0.7) of Tl{U) for non- 
merging galaxies with young GCs (LROO). Also there is a 
factor of ~ 200 difference between the minimum (=0.07) and 




Figure 4. Number of gas-rich ("wet") major mergers that form 
MPCs (solid) and MRCs (dotted) as a function of z. Here wet 
major mergers are those with /g > 0.5 and /m > 0.5. The num- 
bers normalized by their maximum values for < z < 12 are 
shown. 
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Figure 5. The age-metallicity relation for all GCs. The mean 

mctallicity of GCs in age bins arc shown by filled circles. The la 
dispersion in the metallicities of GCs for each age bin is shown 
by an error bar. 



maximum (=15) values of Tl{U) for galaxies with young 
GC candidates in LROO. Therefore 0.005 < a < 0.05 is a 
reasonable range of a in Fm- We try to determine the most 
reasonable values of a and (3 for which both the observed 
net formation efficiency of GCs (McLaughlin 1999) and the 
number fraction of MRCs (Spitler et al. 2007) can be well 
reproduced by the present simulation. 

The net formation efficiency of GCs (e) in the present 
study is defined as follows: 



(7) 



where Mgc and Mstar are total masses of GCs and stars, 
respectively, which are formed in all building blocks in the 
simulation. We define /mrc as follows: 



ATmpc + N^. 



(8) 



where A^mpc and A'mrc arc rmmbers of MPCs and MRCs, 
respectively. Since we find that the model with a = 0.02 and 
/3 = 0.05 can best reproduce the observations, we discuss 
this model in this paper. 
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Figure 6. Metallicity distribution functions (MDFs) for four 
representative massive galaxies with logiQ(^^^) > 13.0: (a) Mb = 

-20.7 mag and B/T = 1.0, (b) Mb = -21.3 mag and B/T = 1.0, 
(c) Mb = -21.3 mag and B/T = 1.0, and (d) Mb = -20.9 mag 
and B/T = 1.0. The mean mctallicitics of MFCs, MRCs, and all 
GCs are shown by solid, dotted, and dashed lines, respectively. 
Although the shapes of MDFs arc different in diflferent galaxies, 
they all clearly show bimodal MDFs. 
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Figure 7. The agc-mctallicity relations (AMRs) of GCs for the 
four galaxies shown in Fig. 6. The mean mctallicitics of GCs are 
shown for the five age bins. In order to show the AMRs more 
clearly, only GCs with ages larger than 10 Gyr arc shown as few 
younger GCs are present. The la dispersion in the metallicities 
of GCs for each age bin is shown by an error bar. 



2.4 Main points of analysis 

The total number of virialized halos at 2: = is 95139, 
among which only 12179 (thus 12.8 %) have GCs that arc 
formed in their building blocks virialized before z = 6 (i.e., 
reionization) . The total number of GC particles at z = 
is 998529 in the present simulation. Fig. 1 shows the large- 
scale structure of GCs that are in the very central regions 
of their host galaxies at z = 0. Although these central GCs 
in galaxies are relatively metal-rich with a mean metallic- 
ity of [Fe/H] = —0.18 for the 12179 galaxies, some fraction 
of them (37.5%) are metal-poor ([Fo/H] < -1). Hero wo fo- 
cus on correlations between physical properties of GCSs and 
those of their host galaxies; we do not discuss the internal 
properties of individual GCSs. 

In investigating GCS properties, we divide GCs into 
MFCs and MRCs according to their metallicities ([Fc/H]): 
those with [Fe/H] < — 1 and with [Fe/H] > — 1 arc defined as 
MFCs and MRCs, respectively. We investigate the physical 
properties of GCSs separately for these two GC subpopu- 
lations. We also discuss the bimodality in the metallicity 
distribution functions (MDFs) of GCSs. Wo show some ex- 
amples of individual representative GCSs in Appendix A. 

Following the morphological classification scheme by 
Simien & de Vaucouleurs (1986), we classify simulated galax- 
ies into three different morphological types: E, SO, and Sp. 
We use the simulated B-band bulgc-to-disk luminosity ratio, 
B/D in the above morphological clEissification. In this paper, 
galaxies with B/D > 1.52, 0.68 < B/D < 1.52, and B/D < 
0.68 axe classified as elliptical (E), lenticular (SO), and spiral 
galaxies (Sp), respectively. 



3 RESULTS 

3.1 Cosmic evolution 

Since our main focus is to discuss the results of correlations 
between physical properties of GCSs and those of their host 
galaxies aX z = 0, we do not discuss the cosmic evolution of 
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Figure 8. The number fractions of GCSs with bimodal MDFs 
(/bimo) a function of Mb for Es (top), SOs (middle), and Sp 
(bottom). The error bar in each bin is based on Poisson a/JV 
statistics. 
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Figure 9. Distributions of tlie simulated GCSs in the /mrc- 
B/T plane (upper) and the /mrc-Me one (lower), where /mrc is 
the number fraction of MRCs in a GCS. Triangles and squares 
represent the observational results by Spitler et al. (2007) and by 
P06, respectively, whereas circles show the average values of /mrc 
in B/T and Mb bins. The la dispersion for each bin is shown by 
an error bar. 
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Figure 10. Distributions of the simulated GCSs on the fy-B/T 
plane (upper) and on the fy-M^ one (lower), where fy is the 
number fraction of young MRCs with ages less than 8 Gyrs. Cir- 
cles show the average values of fy in B/T and Mb bins. The 1<7 
dispersion for each bin is shown by an error bar. 



Table 2. Metallicity-luminosity relations. The observed (o; Peng 
et al. 2006) and simulated (s) coefficients are given for MFC, 
MRC and all GCs, where [Fe/H] = a -t- b Mb. 



MFC 
MRC 
All 



3>0 

-2.771 
-2.452 
-4.854 



bo 
-0.064 
-0.105 
-0.206 



as bs 

-2.414 -0.052 

-2.056 -0.081 

-3.910 -0.160 



GCS. We do however briefly summarize the cosmic evolution 

of GCS properties averaged over all GCs formed in the sim- 
ulation. Fig. 2 shows that the peaks of cosmic GCFRs arc at 
« - 7 for MFCs and at a ~ 4 for MRCs and that the GCFR 
is much higher in MFCs than in MRCs at a > 2, whereas it 
is higher in MRCs for z < 2. Fig. 3 shows that number distri- 
butions of MFCs and MRCs iiorniahzcd by their maximum 
values have two different peaks, which reflect the differences 
in cosmic evolution of GCFRs between MFCs and MRCs. 

Fig. 4 shows that both MFCs and MRCs can be formed 
from gas- rich (/g > 0.5), major (/m > 0.5) mergers for a 
wide range of rcdshifts and that formation rates of MFCs 
and MRCs from wet major mergers peak at 2: ~ 5 for MFCs 
and at « ~ 3 for MRCs. As shown in Fig. 5, the mean 
metalhcities for all GCs (including both MFCs and MRCs) 
steeply increase for ages > 10 Gyr as a result of rapid chem- 
ical enrichment in the building blocks of galaxies. The time 
evolution of the mean mctallicities appears to be much less 
dramatic for ages < 10 Gyr owing to the slower chemi- 
cal enrichment processes in galaxies. The dispersions in the 
mean mctallicities ((7(Fc/H)) arc smaller for larger ages i.e., 
older GCs. We find o-(Fe/H) ~ 1.6 dex for ~ 13 Gyr and 
o-(Fe/H) ~ 0.1 dex for ~ 0.8 Gyr. 

It is interesting to derive the age-metallicity relation of 
GCs with ages t (Gyr) older than 10 Gyr in the present sim- 
ulation at z = in order to compare the simulated relation 
with the observed one when sufficient observational results 
are available. The x-square fit to the simulation data gives 
the following: 



[Fe/H] = 10.86 - 10.831ogiot 



(9) 



The model predicts the presence of young (< 3 Gyr) GCs 
with [Fe/H] > 0, though the number fraction of such GCs is 
quite small. Most of these metal-rich GCs arc found to be 
located in the central regions of galaxies in the simulation. It 
should be however stressed that very old (> 13 Gyr), metal- 
poor ([Fc/H] < —1.6) GCs can also be found in the very 
central regions of galaxies at 2: = 0. 



3.2 Metallicities and ages 

Fig. 6 shows that the GCSs in four early-type galaxies with 
B/T « 1 at « = have bimodal MDFs with peaks below 
and above [Fe/H] = — 1. The locations of the two peaks and 
the number ratios of GCs in each peak arc different between 
different galaxies, which implies that merging and star for- 
mation histories of their building blocks are quite diverse 
in the four galaxies. The peak metallicities of GCSs reflect 
the gaseous metallicities of the building blocks at the epoch 
when strong starbursts were triggered mostly by gas-rich 
("wet"), major merging at high z (> 3). Although these 
four galaxies located in the centers of massive clusters of 
galaxies with Mh > lO^^M© clearly show bimodal MDFs, 
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Figure 11. Distributions of the simulated GCSs in the Mb- 
[Fe/H] plane for all GCs (top), MRCs (middle), and MFCs (bot- 
tom). Here [Fe/H] represents the mean metallicity of a GCS in 
a galaxy. Circles show the average values of [Fe/H] in a Mb bin. 
The 1<T dispersion for each bin is shown by an error bar. For com- 
parison, the observed Mb- [Fe/H] relations of P06 are shown by 
solid lines in the three panels. 
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Figure 12. Distributions of the simulated GCSs in the Sn-Mb 
plane for E (top), SO (middle), and Sp (bottom). Circles show the 
average values of 5n in a Mb bin. Open squares represent the 
observational results by P06, for which we assume that B — V = 
0.9 for all Es to estimate 5n. The Icr dispersion for each bin is 
shown by an error bar. 



less luminous galaxies in the field and groups of galaxies do 
not necessarily show bimodal MDFs in their GCSs. The di- 
verse shapes of MDFs of these galaxies (including those that 
show unimodal MDFs) are briefly discussed in Appendix A. 

Fig. 7 clearly shows that younger GCs arc likely to 
be more metal-rich for all four galaxies and that the age- 
metallicity relations are similar to one another. Since the 
colours of GCs depend on both age and metallicity, these 
results suggest that the colour distributions of the simulated 
GCSs will bo somewhat different from their MDFs. On the 
other hand, the GCSs are dominated by very old GCs in each 
case. We plan to investigate whether the colour distributions 
of GCSs are consistent with observations by combining the 
present GC formation model with a stellar population syn- 
thesis code and the efFects of photometric uncertainty in a 
future paper. 

To investigate the number fractions (/bimo) of galax- 
ies that reveal clear bimodal MDFs we determine /bimo as 
follows: 

/bimo = ,7™° ' (-"^^^ 



where A^bimo and A'gai arc the number of galaxies with obvi- 
ously bimodal GCS MDFs and that of all galaxies, respec- 
tively. In the present study, an MDF with obvious bimodal- 
ity is defined as revealing a clear metallicity peak at both low 
([Fe/H] < -1 for MFCs) and high metallicity ([Fe/H] > -1 
for MRCs). 

Fig. 8 shows that more than 70% of galaxies with 
Mb < — 19 mag have bimodal MDFs in their GCSs regard- 
less of their Hubble types. The number fraction of GCSs 
with bimodal MDFs (/bimo) depends on Mb such that less 
luminous galaxies have smaller /bimo (i.e., less likely to have 
bimodal MDFs). Considering that we do not include GCSs 
with no MFCs in estimating /bimo, the above result means 
that only MFCs arc formed in the building blocks of these 
less luminous galaxies with no bimodal MDFs. Further- 
more gas-rich, major merging between high-density build- 
ing blocks does not occur after their gaseous metallicities 
becomes higher than [Fe/H] ~ — 1 in the formation histories 
of these galaxies. 

Fig. 9 shows that the number fraction of MRCs (/mrc) 
does not depend on B/T. Some very late-type spirals with 
B/T < 0.05 have high values of /mrc (> 0.6) whereas some 
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Figure 13. The same as Fig. 12 but for tlie Vm-Mb plane, wtiere 
Vn is the number of GCs per halo mass. 



ellipticals with B/T fs 1 have low values of /mrc (< 0.1), 
which suggests that the origin of the Hubble types of galax- 
ies are not closely associated with the metal-rich fraction 
(/mrc) in their GCSs. However /mrc docs depend on Mb 
such that /mrc is higher in more luminous gala^xies, though 
the dispersions in /mrc are quite large over a wide range of 
Mb- Galaxies with Mb < —20 mag cover a similar range 
for observed galaxies in the Mb — /mrc plane (Spitler et al. 
2007). The mean value of /mrc for all GCSs in the present 
simulation is 0.36, which is consistent with that found by 
Spitler et al. (2007). We note that the Spitler et al. sample 
includes non-central galaxies. 

The mean formation epochs of MPCs and MRCs for all 
GCs formed in the simulation between 2: = 41 and z = are 
estimated to be 5.7 and 4.3, respectively. These results imply 
typical ages of MPCs and MRCs at a = are 12.7 and 12.3 
Gyr, respectively. This is consistent with the current spec- 
troscopic measurements of extragalactic GCs (e.g. Brodie & 
Stradcr 2006). The mean values of the number fractions of 
young GCs with ages less than 8 Gyrs (/y) is very small 
(<4%) in the present simulation. Fig. 10 shows that fy does 
not depend on B/T. The dependence of fy on Mb appears 
to show no/little trend: although fy appear to peak around 
Mb = —19, this apparent trend is of little statistical sig- 
nificance, owing to the large dispersions. Dispersions in the 
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Figure 14. Distributions of the simulated GCSs in the ^N-Afh 
plane for E (top), SO (middle), and Sp (bottom). Circles with 
error bars show the average values of 5n in Mh bins. The la 
dispersion for each bin is shown by an error bar. 



locations of galaxies in the B/T — fy and the Mb — fy planes 
reflect the diversity in the formation epochs of young MRCs, 
which are mostly via gas-rich major merging. 

Next we investigate correlations between mean metal- 
licities of GCs (simply referred to as [Fe/H]) and Mb of their 
host galaxies. These correlations are called "metallicity- 
luminosity relations" for convenience in the present study. 
The observed metallicity-luminosity relations by P06, after 
transformation from observed colours, can be described as 
follows: 



[Fe/H] = a + 6Mb 



(11) 



The values of the observed coefficients Oo and 60 for MPCs, 
MRCs, and all GCs are summarized in Table 2. The values 
of our simulated Us and 6s using the same relation are listed 
also in Table 2 for comparison. 

Fig. 11 shows that the mean metallicities of all GCs 
correlate well with Mb such that the metallicities are higher 
in more luminous galaxies. The simulated slope of the 
[Fc/H] — Mb correlation for all GCs is consistent with the ob- 
served one by P06 for non-central galaxies. The simulated 
[Fe/H] — Mb correlation for MRCs is also consistent with 
the observed one which suggests that the present formation 
model for MRCs is realistic. The [Fe/H] — Mb correlation for 
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Figure 15. The same as Fig. 14 but for the Vn-Mji plane. 



Table 3. Mean Sn and Vjsi in different Hubble types 



< Sn > 

< > 



Sp 
1.8 
8.3 



SO E 
2.0 4.0 
11.4 11.5 



MPCs is, however, not so consistent with the observed one: 
the slope of the correlation is too flat and the mean metal- 
licitics of MPCs arc systematically higher than the observed 
ones for a wide range of luminosities. One way to solve this 
inconsistency is discussed later in §4.3. 



3.3 Specific Frequencies and 14r 

We investigate specific frequencies (S'n) of GCSs in galaxies 
with different Hubble types. is defined as follows (Harris 
& van den Bergh 1981): 



^^q0.4(M.+16) 



(12) 



where A'gc and Mv arc the total number of globular clus- 
ters in a galaxy and the band absolute magnitude of the 
galaxy, respectively. We also investigate the number of GC 
per unit halo mass V^n, which, following Spitler et al. (2007), 
is defined as follows: 



Fn =iVgc( 



IQIIMb 



(13) 



where Mh is the total mass in a galaxy halo including dark 
matter. The mean values of Sn and Vn in GCSs of galaxies 
of different Hubble types are summarized in Table 3. We 
estimate the la dispersion in St^ and Vn for Mb and Mh bins 
to investigate the statistical significance in the simulated 
correlations of S'n and Vn with Mb and Mh. In estimating 
these dispersions, we use the simulated galaxies with S'n 
{Vn) less than 20 to avoid unreasonably large dispersions 
caused by a very small number of galaxies with unusually 
large Sn (> 50). 

Fig. 12 shows that Sn of E/SOs are typically higher than 
those of spirals (Sp) for a given Mb. Although more lumi- 
nous Es with Mb < —20 have higher Sn than intermediate- 
luminosity Es with — 19 < Mb < —17, there is no such a 
trend in spirals (Sp). Faint galaxies with Mb > — 15 show 
high Sn (> 4) which for E/SO is due to high halo mass- 
to-light-ratios (Mh/La). The mean 5n for Sp, SO, and E 
populations are 1.8, 2.0, and 4.0, respectively. 

Fig. 13 shows that the dependences of Vn on the Hub- 
ble types of galaxies are not as strong as those of S'n: only 
a factor of 2 — 3 difference in Vn between E/SO galaxies and 
spirals for a given luminosity. For a given Hubble type, Vn 
does not depend strongly on Mb. These results mean that 
the numbers of GCs per unit halo mass (Mh) does not de- 
pend strongly on luminosity or the Hubble type of their host 
galaxies. The weak dependence of Vn on Mh was already 
pointed by Bekki et al. (2006, B06), though V^n is estimated 
only for MPCs in B06. 

There is a weak tendency for more luminous Es with 
Mb < —20 mag to have higher Vn, though the dispersion 
in V^ is large. SOs have V^n values > 8 that are significantly 
higher than those of spirals for a wide range of Mb, suggest- 
ing that the higher Vn in SOs is not due to the truncation of 
star formation in spirals as Vn does not change when a disk 
fades via truncation of star formation. These results imply 
that if the observed V^n are typically higher in SOs than in 
spirals, then only a minor fraction of spirals can be trans- 
formed into SOs via truncation of star formation in spirals. 

Fig. 14 shows an "U-shape" distribution in the Sn — Mh 

results for Es, although error bars are quite large (up to 
~ 5 for low-mass halos). This U-shape distribution means 
that Sn is higher in smaller Afh below the threshold halo 
mass (Mh,th) of logj^o-^h.th = lO'^'^M© whereas it is higher 
in larger Mh above Mh,th. This simulated distribution for 
Es appears to be broadly consistent with the latest obser- 
vation by Peng et al. (2008) for GCSs in the Virgo cluster 
of galaxies. This U-shape distribution can be also seen in 
spirals, though it is not so clear in SOs. The Sn — Mh rela- 
tion for halo masses smaller than Mh,th is due mainly to the 
fact that A/h/ie depends strongly on Afh (or on Afe). The 
Sn — Afh relation for halo masses larger than Mh,th is due 
mainly to the fact that numbers of GCs per unit halo mass 
are higher in more massive galaxies. 

As shown in Fig 15, the U-shape distributions are less 
well defined for the Vn — Mh results of galaxies with different 
Hubble types. However, Vn in galaxies with Afh > Mh,th are 
significantly higher in Es, which means higher numbers of 
GCs per unit halo mass in more massive Es. More luminous 
spirals show higher Vn for logMh < lO^^'^M©. It is not so 
clear why SOs have Vn values higher than those of Es and 
Sps at low masses logMh < 1O^°-^M0. 
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Table 4. Mean properties of the simulated GCSs. 

<^f,mpc>^ <2f,mrc>^ < /bimo > ^ < /mrc > £ < fy > ^ < Sj^ > £ < Vn > £ < S^ff > ^ < Sdis > 1 < e > j_ 

5.7 4.3 0.52 0.36 0.04 2.11 11.25 0.84 0.98 0.0013 



The 
The 
c The 
^ The 
e The 
/ The 
f The 
^ The 
* The 
j The 



mean 
mean 
mean 
mean 
mean 
mean 
mean 
mean 
mean 



formation epoch of all MPCs in the simulation, 
formation epoch of all MRCs in the simulation. 

number fraction of galaxies with GCSs showing clear bimodal MDFs. 
number fraction of MRCs in GCSs. 

number fraction of young GCs with ages younger than 8 Gyrs in GCSs. 
Sn (GC number per unit luminosity) in galaxies with GCs. 
Vn (GC number per unit halo mass) in galaxies with GCs. 
ratio of half-number radii of MRCs to those of MPCs in GCSs. 
ratio of velocity dispersions of MRCs to those of MPCs in GCSs. 
net formation efficiency of GCs in galaxies. 
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Figure 16. Dependences of rc,mrc/r-c,mpc (top) re,mrc (middle), 
and re,mpc (bottom) with Mb- The values averaged for GCSs in 
Mb bins are shown by filled circles. The la dispersion for each 
bin is shown by an error bar. 



The results in Figs. 14 and 15 thus show that Sn and 
are fairly high in Es located in the centers of massive 
halos with logMh ~ lO^^'^M©. These results suggest that 
giant elliptical galaxies (and cDs) in the centers of massive 
groups and clusters of galaxies have higher Sn and Vn than 
those in the field and small groups of galaxies. 
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Figure 17. Dependence of ro.mpc (i.e., the half-number radius 
of a GCS) on Mji (halo mass of the host). The values averaged for 
GCSs in Mil bins are shown by filled circles. The Icr dispersion 
for each bin is shown by an error bar. The x-squarc fit to the 
simulation data is shown by a solid line. 



3.4 Dynamical properties 

3.4-1 Structures 

We next investigate the half- number radii of MPCs (r-c.mpc) 
and MRCs (re,mrc) for GCs within the virial radii (rvir) of 
host galaxy dark matter halos: it should be stressed that 
the results below might well depend weakly on the adopted 
initial distributions of GCs within their host halos at high z- 
The spatial distributions of MPCs have valuable information 
of dark matter halos of their host galaxies (B07). In order 
to avoid contributions from intracluster and intragroup GCs 
in groups and clusters of galaxies with rvir > 100 kpc, we 



and clusters. We define Sett to be 



Soff = 



(14) 



We also define the power-law slope (7) and the coefficient 
(Do) of the following relation: 



(15) 



Fig. 16 shows that rc,mpc and re,mrc depend on Mb such 
that both are larger for more luminous galaxies. The mean 
values of re,mpc are larger than 10 kpc for —22 < Mb < —14. 
This is larger than the observed the half-number radius of 
the Galactic GC system which is about 5 kpc (e.g., van den 
Bergh 2000). Although no observational studies have so far 
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Figure 18. Dependences of <Tmrc/o"mpc (top) (Jmrc (middle), and 
o"mpc (bottom) on Mb in the simulation. The values averaged for 
GCSs at each of the six Mb bins are shown by triangles. The la 
dispersion for each bin is shown by an error bar. 



which means re.mpc = 0.20(^^)°'^^ kpc (i.e., 7=0.18). The 
following equation would be more useful for making an esti- 
mation of Mh for a galaxy by measuring re,mrc of the GCS: 



:4.5-^5.21ogio(- 



kpc 



(17) 



Since the present model based on a dissipationless sim- 



ulation can overestimate re 



it might be better to use ob- 



servations in order to determine the zero-point of the above 
power-law relation. If we use observations of the Galactic 
GCS (e.g., van den Bergh 2000 and references therein) and 
the total mass of the Galaxy (Wilkinson & Evans 1999), 



then the power-law relation is 



kpc or Mh 



5-0(2x1012Mo 



2 X 10^^(- 



5kpc 



M©. The derived Mh from 

-Mh relations can be compared with Mh derived from 
kinematics of GCSs and halo field stars (e.g., Romanowsky 
2006). 



3.4-2 Kinematics 

A growing number of observational data sets on the kino- 
matical properties of GCSs have been recently accunmlatcd 
(e.g., Romanowsky 2006). We investigate correlations be- 
tween velocity dispersions (a) for MPCs ((Tmpc) and MRCs 
(cTmrc) and Mb of their host galaxies. We first estimate (Tmpc 
and (Tnirc of a GCS for each of the three projections {x-y, x- 
z, and y-z) by using line-of-sight-velocities of all GCs within 
the virial radius of the halo. We then malce an average for 
the three projections and determine one-dimensional veloc- 
ity dispersions of CTmpc and Cnirc- We also investigate corre- 
lations between Mb and Sdis, where Sdis is defined as: 



(Tmrc 



(18) 



investigated the correlation of Te.mpc with Mb, the above 
results imply that the simulation significantly overestimates 
ro,mpc: this may be true for MRCs as well. In the present 
model, all GCs within rvir (virial radius) of the dark matter 
halo of their host galaxy are used for estimation of re.mpc- 
This way of estimating re,mpc would contribute significantly 
to the possible overestimation of re.mpc of GCSs in the sim- 
ulation. 

We find Seff is significantly less than 1 for a wide range 
of luminosities, which means that the distributions of MRCs 

arc more compact than those of MPCs for most GCSs. The 
mean Seff is 0.84 in the present model, which can not be 
currently compared with observations owing to the lack of 
observational studies of Seff in galaxies. The smaller mean 
value of Seff and the weak dependence of Sefi on Mb can be 
tested against future observational studies. We also suggest 
that the slopes of rc,mpc — Mb and rc,mrc — Mb relations 
can be used to constrain theoretical models of GC formar- 
tion based on hierarchical clustering scenarios, because these 
depend on merging histories of galactic building blocks that 
form GCs. 

Fig. 17 shows re,mpc of GCSs correlate with their host 
halo mass such that re.mpc is larger for larger Mh. The x- 
square fit to the simulation data is: 



are higher in 

more luminous galaxies. The essential reason for this depen- 
dence is described as follows. Since the GCSs investigated 
in the present study are for galaxies located in the centers 
of dark matter halos (i.e., GCSs in satellite galaxies are not 
investigated), more luminous galaxies are more likely to be 
embedded in more massive halos (i.e., larger Mh). GCs fol- 
low structures and kinematics of underlying dark matter ha- 
los so that the mean velocity dispersions of GCSs are deter- 
mined mainly by masses of their halos (and by their spatial 
distributions). Therefore the velocity dispersions of GCSs 
are higher in more luminous galaxies. The difference in Sdis 
for different Afe is quite small and the mean value of Sdis is 
0.98. These results imply that velocity dispersions are not 
so different between MPCs and MRCs in galaxies. 

Owing to the coUisionless nature of the present simu- 
lation, the ratios of maximum rotational velocities to cen- 
tral velocity dispersions (V^/f) in GCSs appear to be low 
iy jo < 0.3) for most galaxies. These results are in a striking 
contrast with our previous simulations (Bekki et al. 2005) in 
which some GCSs in the remnants of disk-disk major merg- 
ers show largo V/a (> 0.5). Owing to the lack of extensive 
statistical studies on V/ a of GCSs in galaxies, it is not clear 
whether the above results are generally consistent with ob- 
servations or not. 



Formation of GC systems 13 



N g _ 




10 



-14 



-18 

Mb (mag) 



-20 



-22 



Figure 19. Distributions of the simulated galaxies with GCs in 



B/T plane (upper) and Zy 



Mb plane (lower) 



the ^vifjgc 

where Zvir.gc is the mean z^ir of building blocks hosting GCs in 
a galaxy. Circles show mean values in B/T and Mb bins. The 
la dispersion for each bin is shown by an error bar. No galaxies 
have Zvir.gc < 6, because formation of GCs in the present model 
is truncated for galaxies with z^ir < 6. 



3.5 Mean properties 

Lastly, we briefly summarize the key physical properties of 
GCs averaged over all GCSs in the simulation in Table 4. 
The columns are the mean properties of a;f,mpc (column 1), 

Zf.mrc (2), /bimo (3), /n.rc (4), /y(5), Sn (6),' Vn (7), Seff (8), 

Sdis (9), and e (10). In estimating these mean values, we do 
not include galaxies with no GCs. If wc include galaxies with 
no GCs, the mean values would be significantly changed, in 
particular, for and Vn. 



4 DISCUSSION 

4.1 Bimodal colour distributions in GCSs 

The origin of the observed bimodal colour distributions of 

GCSs in elliptical galaxies have long boon discussed in the 
context of different formation scenarios of elliptical galax- 
ies, such as two-phase collapse at high redshift (Forbes et 
al. 1997), accretion and stripping of low-mass galaxies with 
GCs (Cote et al. 1998), and gas-rich major mergers (Ashman 
& Zopf 1992). B02 was the first to investigate the bimodal 
colour distributions in a quantitative manner based on the 
results of a semi-analytic model of galaxy formation. B02 
showed that the observed bimodality can be reproduced, if 
MRCs are formed during dissipative merger events at high 
redshift and that the formation of MFCs arc truncated at 
« ~ 5. B02 however did not constrain the truncation mech- 
anism. 

Since age differences between MFCs and MRCs are 
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Figure 20. Distributions of the simulated galaxies with GCs in 
the Mjj/Lb — Mb plane, where Lb is the total B— band luminosity 
of a galaxy, for E (top), SO (middle), and Sp (bottom). Circles 
show mean values in Mb bins. The la dispersion for each bin is 
shown by an error bar. 



rather small (< 1 Gyr) in the present simulation, we can 
discuss the simulated MDFs of GCSs in terms of the colour 
bimodality of GCSs. In the present model, all GCs (i.e., 
both MFCs and MRCs) in a galaxy originate from low-mass 
building blocks at z > 3, whether they are isolated or in 
merging galaxies. Therefore the origin of MRCs in a galaxy 
is not necessarily associated with the past major merger 
events of a galaxy. The peak metallicity in the MDF for 
MFCs (MRCs) in a galaxy reflects the highest GCFR in 
the galaxy's building blocks that have stellar metallicities 
lower (higher) than —1 at high z. The mean metallicity 
of MFCs (MRCs) in a galaxy is determined by the mean 
stellar metallicity of the more metal-poor building blocks 
with [Fe/H] < -1 ([Fe/H] > -1). Although MRCs are 
formed later in more metal-rich building blocks of galaxies, 
there is only a slight difi'erence in formation redshifts be- 
tween MFCs and MRCs. The major merger events in which 
MRCs are formed are typically z > 3 whereas the last major 
merger events which determine the final morphological types 
of galaxies can happen later at lower redshifts. Therefore the 
last major merger events responsible for elliptical galaxy for- 
mation are not necessarily associated with the formation of 
MRCs , which is in contrast to the scenario proposed by Ash- 
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man & Zepf (1992) but consistent with the so-called 'damp' 
merger interpretation of Forbes et al. (2007). 

In order to show the key requirements for the formation 
of bimodal MDFs in GCSs, we have run two comparative 
models with a = 0.1 and /3 = 0.05 and a = 0.02 and (3 = 
0.25 (see section 2.3). For these models, the dependences of 
the GCFR on fm and /g are assumed to be very weak. We 
find that GCSs in these models show much reduced MDF 
bimodality in comparison with the present standard model 
with a = 0.02 and /3 = 0.05 shown in Table 1 (the results 
of the models are shown in the Appendix B). These results 
clearly show that a strong enhancement in the GCFR during 
violent merging (in particular, major merging) and gas-rich 
phases (i.e., at high redshifts) is important for the formation 
of the bimodality. 

The origin of the metal-poor GC peaks is due mainly to 
the strong enhancement of the GCFR in gas-rich and high 
gas-density building blocks at high-a; whereas that of the 
metal-rich ones is due mainly to the strong enhancement of 
the GCFR in more metal-rich and gas-rich building blocks 
that experience violent merging at somewhat later epochs. 
The truncation of GC formation via reionization is neces- 
sary for the present model in order to prevent too many 
GCs from forming in low-mass galaxies. The epoch of reion- 
ization, however, appears not to determine the presence of 
bimodality itself. Although differences in the merging his- 
tories of halos between different galaxies can introduce a 
diversity in GCS properties, such differences appear to have 
no direct link with the resulting bimodal MDFs of the GCSs. 

As shown in Fig. 8, not all of galaxies have GCSs with 
bimodal MDFs. In order to understand why these galaxies 
have GCSs without bimodality, wc have investigated forma- 
tion histories of GCs in these galaxies (the results for some 
galaxies are shown in the Appendix B). We find that these 
galaxies have only either metal-poor or metal-rich peaks in 
their MDFs owing to single burst events of GC formation. 
Galaxies with no, or little, enhancement in the GCFR for 
metal-poor GCs and strong enhancement of the GCFR for 
metal-rich GCs have no clear bimodality in their MDFs. 
On the other hand, those with strong enhancement of the 
GCFR for metal-poor GCs at high z and no enhancement of 
the GCFR for metal-rich GCs do not show the bimodality 
either. The efficient formation of MFCs is highly unlikely for 
galaxies with lower «vir whereas MRCs are highly unlikely to 
be formed in gas-poor merging at later epochs. The galaxies 
having GCSs without bimodal MDFs thus have either lower 
ZviT or few events of gas-rich merging and accretion of more 
metal-rich galactic building blocks. 

The present model has shown that a significant fraction 
(/bimo ~ 0.2) of very late-type spirals with small or no bulges 
{B/T < 0.05) have bimodal MDFs in their GCSs. This is 
because more metal-rich ([Fe/H] > —1) building blocks of 
late-typo spiral galaxies can form GCs at high z which are 
later accreted during mergers. Therefore the origin of MRCs 
(and thus bimodal MDFs) in late-type spiral galaxies are not 
necessarily associated with the formation of galactic bulges 
via early major merger events in the present model (see B02 
for a brief discussion on this issue). If most bulge- less spirals 
are found to contain very few MRCs (as appears to be the 
case for M33; Chandar et al. 2006) then the model proposed 
here may require significant alteration. 



4.2 Specific Frequencies 5n and Vn 

Previous observations have revealed that the Sn of GCSs de- 
ponds on the Hubble type and luminosity of the host galaxy 
(e.g., Harris 1991; but see also Spitler et al. 2007 for an al- 
ternative view). Although the observed trends of Sn with 
galactic properties (e.g., the bimodality in the A/v — Sn 
relation) are well reproduced by BOG, other key observa- 
tions of Sn have not been discussed. In the present study, 
we have shown that (i) 5n of GCSs in luminous galaxies 
with Mb < — 19 mag arc significantly higher in early-type 
(E/SO) galaxies than late-typo (Sp) ones, (ii) SOs typically 
have higher than spirals, and (iii) low-luminous galax- 
ies with Mb > —15 mag have higher Sn, regardless of their 
Hubble types. These results are qualitatively consistent with 
previous observations by Harris (1991), Forbes (2005), and 
Aragon-Salamanca et al. (2006). So far wc have not dis- 
cussed the physical reasons for the above dependences of Sn 
on galactic properties. 

Fig. 19 shows that the mean values (zvir.gc) of viri- 
alization redshifts (zvir) of galactic building blocks that 
form GCs arc significantly higher in early-type galaxies 
with B/T > 0.8 than in late-type ones. This means that 
GCFRs in building blocks of early-type galaxies are higher 
owing to their higher mass- densities with higher z^ii. Be- 
cause GCFRs depends strongly on z^ir through the term of 
Fz (see equation 6) the higher «vir,gc in early-type galaxies 
is one reason for the observed higher Sn of their GCSs. Fur- 
thermore Fig. 19 shows that Zvir.gc is higher in more lumi- 
nous galaxies, which can explain why more luminous early- 
type galaxies with Mb < —19 mag are more likely to have 
higher Sn in comparison with less luminous ones with —19.0 
< Mb < -17 mag seen in Fig. 12. 

Fig. 19 shows that low-luminosity galaxies do not have 
higher ^vir.gc, which implies that higher Sn in low-luminosity 
galaxies seen in Fig. 9 can not be understood in terms of 
2vir,gc. B06 suggested that the origin of the observed higher 
Sn in low-luminosity galaxies is due to higher mass-to-light- 
ratios {M/L ~ Mh/I/B) of these galaxies. Fig. 20 shows 
that Mh/I/B steeply depends on luminosity in the sense that 
Mh/Ls is higher in less luminous galaxies. Given that GC 
numbers per unit halo mass (Vjsi) in galaxies only weakly de- 
pend on the halo masses of the galaxy (B06; see also Fig. 13 
in the present paper), the above result implies that the 
higher Sn in low-luminosity galaxies are due mainly to their 
higher Mh/Lb ratios owing to the relation Sn = ^n(^). 

4.3 Metallicity-luminosity relations 

Although B02 predicted that the mean colours of both 
MFCs and MRCs only weakly correlate with the total lumi- 
nosities of their host galaxies, these predictions are incon- 
sistent with latest observations by Strader et al. (2004) and 
P06. The present model has successfully reproduced reason- 
ably well, the observed [Fe/H] — Mb relations for MRCs 
and all GCs. The derived slope in the [Fe/H] — Mb relation 
for MFCs in the present model is, however, not consistent 
with the observed one. The differences between B02 and the 
present study are caused by the differences in the adopted 
models of GC formation between these two studies. Thus 
more sophisticated models of GC formation in galaxies at 
high z need to be developed so that all three [Fe/H] — Mb 
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Table 5. Qualitative assessment of model successes and failures. 



Items 


Consistency ^ 


Comments 


Ages 


o 


Data indicate very old ages for both MRCs and MPCs. 


Bimodal MDFs 


o 


Requires predicted MDFs in observed colours. 


/bimo — A^B 


o 


Broadly consistent but data very limited. 


/mrc — Mb 


o 


Consistent with data for luminous galaxies. 


fy 


o 


Broadly consistent but data very limited. 


Mb - [Fe/H] 


o 


Less consistent for MPCs. 


Sn - B/T 


o 


Consistent in terms of higher Sj-j in Es. 


Sn - Mb 


A 


U-shape of data not clearly seen. 




A 


Data show a flatter slope. 




A 


Data for the Galaxy is much smaller. 



" O (^) means that simulations are broadly consistent (inconsistent) with observations. A means that simulations are only partly consistent 
with observations. 



relations (MPCs, MRCs, and all GCs) can successfully re- 
produce the observations. 

B07 showed that if Ztmn ~ 10, the observed [Fe/H] -Mb 
relation for MPCs is better reproduced. However the simu- 
lated e (i.e., the formation efficiency of GCs) in the model 
with ztrun = 10 in the present study is too small (0.00077) to 
be consistent with the observed one by McLaughlin (1999). 
The present model seems to have difficulties in explaining 
self-consistently both the observed [Fe/H] — Mb relation for 
MPCs and e. One possible way to solve this problem would 
be to investigate models with higher ztmn (=6 — 10) in which 
Ceff depends more strongly on z than in the present model 
so enough GCs can be formed at z = 6 — 10. It is, how- 
ever, unclear whether the models with higher Ztmn and C'ctt 
can reproduce other key observations such as the 5n — Mb 
relation. 



4.4 Success and failures of the present model 

Although we have so far presented variously different physi- 
cal properties of GCSs in galaxies, only some of them can be 
directly compared with observations: for example, the simu- 
lated /mrc — Mb relation can be compared with observations 
by P06 whereas the simulated Sdis — Mb one can not owing 
to the lack of observational data. It is, however, important to 
check whether the simulated properties are consistent with 
the corresponding observations. In Table 5 we summarise 
what we regard as the relative successes and failures of the 
current model in a qualitative sense. For any given physi- 
cal property there are subtle, and sometimes large, differ- 
ences between the observational data and the corresponding 
property predicted by the model. The references of these 
observational properties are already given in the Introduc- 
tion section §1 (e.g., Brodie & Strader 2006) and some of 
the observational data are presented within figures of this 
paper. Less consistent results (e.g., the Mb — [Fe/H] relation 
for MPCs) imply that a more sophisticated model for GC 
formation would be required for more successful modeling. 



5 CONCLUSIONS 

We have investigated the structural, kinematical, and chem- 
ical properties of globular cluster systems (GCSs) in galaxies 



in a self-consistent manner based on high-resolution cosmo- 
logical N-body simulations combined with a semi-analytic 
model of galaxy and globular cluster (GC) formation. We 
have adopted a number of assumptions on formation efficien- 
cies of GCs which depend on the physical properties of their 
host galaxies (e.g., gas mass fraction). We investigated cor- 
relations between physical properties of GCSs and those of 
their host galaxies for ~ 10^ simulated central halo galaxies 
for MPCs with [Fe/H] < -1 and MRCs with [Fe/H] > -1. 

We find: (1) The majority (~ 90%) of GCs currently in 
halos of galaxies are formed in low-mass galaxies ai z > 3 
with mean formation redshifts of MPCs and MRCs being 5.7 
and 4.3, respectively. This corresponds to 12.7 and 12.3 Gyrs 
in lookback time. The majority of MPCs are formed in low- 
mass galaxies that are virialized well before reionization {z > 
6) and thus have higher mass densities. MRCs are formed 
slightly later not only within high z major mergers between 
high-density galaxies but also within gas-rich isolated ones. 

(2) About 50 % of galaxies with GCs show clear bi- 
modalities in their MDFs, though less luminous galaxies 
with Mb fainter than —17 are much less likely to show the 
bimodalities owing to no or few MRCs. The age differences 
between MPCs and MRCs are quite small (< 1 Gyr) in most 
galaxies. The origin of the simulated bimodality in MDFs is 
due to strong dependences of GCFE on /g and /m (i.e., 
higher GCFE in more gas-rich high-z galaxies and in merg- 
ers with larger mass ratios of the merging two galaxies). 

(3) The number fractions of MRCs (/mrc) range from 
to almost 1 with an average of 0.4. The /mrc in galaxies does 
not depend on their Hubble type (i.e., bulge-to-disk-ratios), 
which implies that the formation of MRCs is not necessarily 
associated with bulge formation. The /mrc are likely to be 
smaller for less luminous galaxies. 

(4) The Sn of GCSs are typically higher in ellipticals 
than in spirals, and in low-luminosity galaxies with Mb > 
— 15 regardless of their Hubble types. The mean Sn for Sp, 
SO, and E populations are 1.8, 2.0, and 4.0, respectively. 

(5) The number of GCs per halo mass (Vn) does not 
depend as strongly on the luminosity or the Hubble type 
of the host galaxy as Sn does, which suggests that the GC 
number per unit mass is similar between different galaxies. 
Vm is, however, likely to be higher in luminous ellipticals 
with Mb < -20. 

(6) Although there are no significant differences in Vn 
between spirals and SOs for luminous galaxies (Mb < —20), 
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SOs are more likely to have higher Vn than spirals for less 
luminous galaxies with Mb > —18. These results suggest 
that only luminous SOs with moderately high Sn can be 
transformed spirals via truncation of star formation and the 
resultant disk fading, because Vn does not change during 
disk fading. 

(7) The mean metallicities of GCs ([Fe/H]) for MFCs 
and MRCs depend on Mb of their host galaxies such that 
they are higher in more luminous galaxies, though the de- 
pendence for MFCs is weak. Although the observed corre- 
lation of Mb — [Fe/H] for MRCs can be well reproduced by 
the present model, that for MFCs is not consistent with the 
simulated one. The observed correlation between the mean 
metallicities of all GCs and Mb can be well reproduced by 
the present model. 

(8) Spatial distributions of MRCs are more compact 
than those of MFCs with rc.mrc ~ 0.84re,mpc • The re,mpc 
depends strongly on halo mass Mh such that rc,mpc oc Afh"'^* 
(or Mh oc re^'^) which implies that re,mpc can be used for 
estimating the total masses of dark matter halos. 

(9) There is no significant difference in velocity dis- 
persion between MFCs and MRCs. Velocity dispersions are 
larger in more luminous galaxies both for MFCs and MRCs 
regardless of their Hubble types. 

(10) The physical properties of GCSs such as MDFs, 
/mrc, Sn, and Vn are quite diverse between different galax- 
ies, depending on the virialization redshifts of their building 
blocks and subsequent merging, star formation, and chemi- 
cal evolution histories of the building blocks. 
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Figure Al. SFR (dotted) and GCFR (solid) as a function of 
z in an early-type galaxy (Gl) located in the central region of 

a massive cluster of galaxies with = 6.3 X 10"'^''Mq. Gl has 
a B/T = 1.0, Mb = -20.7 mag, and My = -21.5 mag. For 
comparison, the SFR and GCFR normalized by their maximum 
values are shown. 
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Figure A2. Solid and dotted lines represent the MDF of GCs 
in the galaxy Gl and that of all field stars (FS) that were formed 
in galaxies previously hosting the GCs of Gl. Note that only the 
GCS shows a clear bimodal MDF. 
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APPENDIX A: INDIVIDUAL CASES 

Here we briefly discuss some results of GCSs in individual 
model galaxies (Gl and G2). Fig. Al shows an elliptical 
galaxy (Gl) at the center of a massive cluster of galaxies 
with Afh = 6.4 X IO^'^Mq which has two strong peaks in the 
GCFR at high redshifts, z ^ 7 and ~ 10. On the other hand, 
SFRs of building blocks of Gl (including Gl) show multiple 
peaks at 2 < z < 6, which means that formation of field 
stars can continue to be quite active till relatively recently. 
For Gl, both the vast majority of field stars and GCs are 
formed till z = 2. Fig. A2 shows that although the MDF of 
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Figure A3. The same as Fig. Al but for a spiral galaxy G2 with 
Mh = 4.1 X 1O"M0. G2 has a B/T = 0.13, Mb = -18.9 mag, 
and My = —19.6 mag. 
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Figure A4. The same as A2 but for a spiral galaxy G2. 

GCs shows clearly a bimodality, that of field stars does not. 
Given that some fraction of these field stars located m the 
outer part of Gl can finally become halo field stars in Gl 
in the present model, the above result suggests that stellar 
halos of ellipticals can have unimodal MDFs as opposed to 
bimodal ones seen in their GCSs. 

Fig. A3 shows that in a simulated late- type spiral (G2), 
the vast majority of GCs are formed before 2 = 3 whereas 
field stars can continue to form till z = 0, which is in sig- 
nificant contrast with the results for the cluster ellipticals 
shown in Fig. Al. It should be stressed here that although 
building blocks of G2 are virialized before z = 6, inten- 
sive GC formation happens around ~ 3. This means that 
GC formation in G2 is triggered by merging of the build- 
ing blocks (with Zvir > 6) at relatively later redshifts. As 
shown in Fig. A4, G2 also shows a lower peak metallicity 
for MRCs ([Fe/H] ~ -0.6) and a smaller number of MRCs. 
The peak metallicity for metal-poor field stars is almost co- 
incident with that for MPCs, which suggests that the MDFs 
are not so different between MPCs and halo field stars in G2. 
It should be stressed here that some galaxies with Mb, Mh, 
and Hubble types similar to G2 have shapes of MDFs in 
GCSs significantly different from that of G2. The variety in 
merging histories between different galaxies can cause these 
differences in MDFs. 

Fig. A5 shows that the MDFs in spirals embedded in 
dark matter halos with different halos masses have variously 
different MDFs: spirals labeled as (a), (c), and (j) show re- 
markably strong peaks for MPCs owing to a larger lumber 
of MPCs whereas those labeled as (b), (e), (h), and (i) show 
clearly multiple (triple) peaks in MDFs. Spirals labeled as 
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Figure A5. A collection of MDFs of 16 GCSs in spiral galaxies 
with B/D < 0.68 and with different Afjj. The 16 spiral galax- 
ies are selected from those embedded in dark matter halos with 
11.5 < logio(^) < 12.0. Here GCSs are selected such that 
just one CCS can be allocated in each of 16 equally spaced mass 
bins in logj^oMj^ over the above mass range. GCSs in lower Mh 
arc given earlier characters in alphabetical orders; the GCSs in 
the lowest and highest Mj^ are labeled as (a) and (1), respec- 
tively. The threshold mctallicitics between MFCs and MRCs (i.e., 
[Fe/H = —1) are shown by dotted lines for comparison. 
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Figure A6. The same as Fig. A5 but for 16 GCSs in carly-typc 
(E and SO) galaxies with B/D> 0.68 and with different Mh. The 
16 early- type galaxies are selected from those embedded in dark 
matter halos with 11.5 < logio(^) < 13.0. 
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Figure Bl. The same as Fig. A2 but for the model with a = 0.1 
and 13 = 0.05. 
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Figure B2. The same as Fig. A2 but for the model with a = 0.02 
and P = 0.25. 



(f) and (1) can be considered to be showing clearly bimodal 
MDFs in their GCSs. The relative contributions of MRCs to 
MDFs arc clearly different between different GCSs in these 
16 spirals. Fig. A6 shows that MDFs in ellipticals embed- 
ded in dark matter halos with different halos masses have 
variously different MDFs: a larger number of GCSs appear 
to have bimodal MDFs in ellipticals than in spirals. These 
diversities in the shapes of MDFs in galaxies with differ- 
ent luminosities and Hubble types reflect the diversities in 
merging and star formation histories between the galaxies. 



APPENDIX B: COMPARATIVE MODELS 

We ran models with different a and f3 in order to more 
clearly understand the origin of the bimodal MDFs seen in 
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Figure B3. The same as Fig. Al but for galaxy G3. 



Formation of GC systems 19 







GC 























[Fe/H] 

Figure B4. The same as Fig. A2 but for galaxy G3. 



owing to the too many metal-rich GCs. Other galaxies hav- 
ing GCSs without the bimodality (which are mostly less- 
luminous systems) do not have metal-rich peaks owing to 
the lack of strong enhancement of GCFR later in their evo- 
lution. The galaxy labeled as (j) with no strong metal-rich 
peak in Fig. A6 is a good example of such galaxies. 



the present simulation. Fig. Bl shows the MDF of the galaxy 
Gl (shown in Figs. Al and A2) for the model with a = 0.1 
(i.e, five times larger than that used in the model described 
in the main text) and (3 = 0.05. For this model with very 
weak dependence of GCFR on /m owing to a larger value of 
a, strong enhancement in GCFR during major merging can 
not happen. The MDF, which can be compared with that 
in Fig. Al, does not clearly show the metal-rich peak. This 
means that strong enhancement in the GCFR during violent 
merging is important for the appearance of bimodality. 

Fig. B2 shows the MDF of Gl for the model with 
a = 0.02 and /3 = 0.25 (i.e., five times larger than that 
used in the model described in the main text). For this 
model with very weak dependence of GCFR on /g, strong 
enhancement in GCFR during gas-rich phases (thus mostly 
at high-redshifts) can not happen. Although the metal-poor 
peak can bo barely seen, the bimodal MDF in this model 
becomes much less significant than that in the model shown 
in Fig. Al. This result suggests that strong enhancement 
of the GCFR in gas-rich phases of galaxies (thus at high- 
redshifts) can be also important for the formation of the 
bimodality. Since we assume that galaxies with higher Zvir 
have higher surface mass densities, stronger enhancement 
of GCFR in more gas-rich galaxies at higher z means high 
GCFR in galaxies with higher surface densities of gas. These 
results thus suggest that the origin of the bimodality can 
result from stronger enhancement of the GCFR in more vi- 
olent merging and and more gas-rich (and higher surface gas 
density) galactic building blocks. 

As shown in Fig. 8 and A6, not all of the simulated 
GCSs have bimodal MDFs. It is therefore important to 
clarify the reasons why these GCSs do not have bimodal 
MDFs. By investigating both MDFs and star formation his- 
tories (SFRs) in galaxies having GCSs without bimodality, 
we find that most of GCSs without bimodality either show 
only metal-poor peaks or show only metal-rich ones owing to 
very strong enhancement of the GCFR during single burst 
epochs of GC formation. Fig. B3 and B4 show the SFH and 
the MDF, respectively, for the galaxy G3 with B /T = 0.7 
and Mh = 2 x 10^^ Mq. Owing to the very strong burst of 
GC formation during late major merging around z — 1 for 
this galaxy G3, too many metal-rich GCs can be formed. 
This galaxy does not experience strong enhancement of the 
GCFR before z — 6 so that the metal-poor peak in the 
MDF is not clear. As a result of these, G3 does not show 
the bimodality in the GCS. G3 is the representative case for 
galaxies which do not show the bimodality in their GCSs 



